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Abstract 

F ^ 1 A new nonlinear Schrodinger equation is obtained explicitly from the 

fractal Brownian motion of a massive particle with a complex-valued dif- 
fusion constant. Real-valued energy (momentum) plane wave and soliton 
' solutions are found in the free particle case. The hydro-dynamical model 

analog yields another (new) nonlinear QM wave equation with physically 
meaningful soliton solutions. One remarkable feature of this nonlinear 
Schrodinger equation based on a fractal Brownian motion model, over all 
, the other nonlinear QM models, is that the quantum-mechanical energy 

■ functional coincides with the field theory one. 

(N 
O 

(N ■ 1 Introduction 

O 

r- 1 ■ The theoretical study of quantum chaos has been developed mainly in two ar- 

eas: The phenomenological characterization of the spacing of the energy levels 
of bound and quasi-bound quantum physical systems, whose main analytical 
tool is the random matrix theory ||, and the semi-classical limit of chaotic 
classical systems |?| . The semi-classical approach pretends to seek solutions of 
the Schrodinger equation and to read in the wave functions any fingerprints 
of classical chaos. Due to the linearity of the Schrodinger equation there is no 
place where the sensibility to the initial conditions can be made manifest, which 
is present in nonlinear chaotic systems. The Riemann zeta function has been 
considered as a unifying link between those two approaches ||. 

We believe that quantum chaos is truly a new paradigm in physics associated 
with non-unitary and nonlinear QM processes based on non-Hcrmitian operators 
(implementing time symmetry breaking). This chaotic behavior stems directly 
from the nonlinear Schrodinger equation without any reference to the nonlinear 
behavior of the classical limit. See ||. For this reason, the genuine quantum 
chaos should be exhibited only by systems whose behavior is correctly described 
by a nonlinear Schrodinger equation. 

The nonlinear QM has a practical importance in different fields, like con- 
densed matter, quantum optics and atomic and molecular physics; even quan- 
tum gravity may involve nonlinear QM. Another important example is in the 
modern field of quantum computing. If quantum states exhibit small nonlin- 
earities during their temporal evolution, then quantum computers can be used 
to solve NP-complete (non polynomial) and #P problems in polynomial time. 
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Abrams and Lloyd [g6| proposed logical gates based on non linear Schrodinger 
equations and suggested that a further step in quantum computing consists 
in finding physical systems whose evolution is amenable to be described by a 
NLSE. 

On other hand, we consider that Nottale and Ord's formulation of quantum 
mechanics Q from first principles based on the combination of scale relativity 
and fractal space-time is a very promising field of future research. In this work 
we extend Nottale and Ord's ideas to derive the nonlinear Schrodinger equation. 
This could shed some light on the physical systems which could be appropriately 
described by the nonlinear Schrodinger equation derived in what follows. 

The contents of this work are the following. In section || we derive different 
nonlinear Schrodingcr-like equations starting from purely hydro-dynamical con- 
siderations. In section [| a review of the derivation of the Schrodinger equation, 
based on Nottale and Ord's [[[) model of QM as a fractal Brownian motion of 
a particle zigzagging back and forth in space-time, is presented. In section || 
we derive the nonlinear Schrodinger equation from an extension of the Nottale's 
approach to the case of a fractal Brownian motion with a complex diffusion 
constant. In section ^ real- valued energy solutions of the nonlinear Schrodinger 
equation are proposed. In the final section ^|, we summarize our conclusions and 
include some additional comments. 



2 Nonlinear Srhodinger equations based on hy- 
drodynamics 

In this section we will write down two NLSE (nonlinear Srhodinger equations) 
using the hydro-dynamical models of QM j|, [l0| . The first equation is based by 
adding a hydrostatic pressure term to the Euler-Lagrange equations Q and the 
second equation is obtained by adding, instead, a kinematic pressure term. As 
far as we know, this second equation has not appeared in the literature before. 

The hydrostatic pressure experienced by a fluid element at a point r due to 
the force of gravity is given by the Euler equation — Vp = pg; V is the ordinary 
gradient, p the pressure, p the density and g the acceleration of gravity. For 
example, if the density and acceleration are uniform one can integrate such 
equation and arrived at p = pgx giving the pressure at a given depth x. 

The author Q proposed to establish the QM analog of the Euler equation 
by relating the density p to the quantum mechanical probability density ip*ip 
and by integrating the equation. Setting p — b a mass-energy parameter 

and the particular case that p = p. then one has that the hydrostatic potential 
is given by the integral 

b [ g(x) -dr = -b I — ■ dr = -bin — = -bln(ib*ib), (1) 
J J P Po 

setting pq = 1. This is the nonlinear potential energy induced from a hydrostatic 
pressure term. 
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It is important to normalize the logarithms by a constant which we set to 
unity for convention. — bhx{tj)* ij)) has energy units. From now these logarithmic 
terms are normalized that way. 

The hydrostatic pressure term [Q has energy units and explains in a straight- 
forward fashion the nonlinear term (nonlinear potential) added to the standard 
Schrodinger equation by Bialynicki-Birula and Mycielski || long ago. The pa- 
rameter b has units of mass (energy), so the nonlinear wave equation is given 
by H after adding the nonlinear potential — bln(ip*ip). 

The Birula-Mycielski NLSE for a particle is 

Bib ft 2 

th-^^~—V 2 i; + U^-b[Hr^. (2) 

A derivation of this equation from the Nelson stochastic QM was given by 
Lemos (|L6| p. 615 and ]l7|). As interesting as this equation may be there are 
some problems. Such equation does not obey the homogeneity condition (see 
Weinberg |Q) which says that if \ip) represents a physical state, the rays \Xip} 
must also represent the same physical state, for any complex constant A. But 
equation (||) is not invariant under ip — > Xib because the logarithmic nonlinear 
potential breaks such homogeneity, this NLSE is not scaled by A. 

Another problem is that plane wave solutions to equation (^|) do not seem to 
have a physical interpretation due to extraneous dispersion relations. Only the 
soliton solutions were physically meaningful Upper limits on the values of 
the parameter b had been found to be |l9j b < 3 • lCU 15 eU, which correspond to 
an electron soliton width of 3 mm The smallness of b is itself no reason to 
disregard equation (^) as physically relevant. For the authors, another problem 
with the logarithmic nonlinear potential term is that the hydrostatic pressure 
term in the NLSE is given by an explicit function of both ip and it complex 
conjugate ip*. It is desirable to write a NLSE solely in terms of the ib variable, 
or ip* , but not combined. Thus, another new NLSE can immediately be writ- 
ten and/or modified by adding the kinetic-pressure terms to the Euler-Newton 
hydro-dynamical equations of motion, i.e. by adding the term (l/2)pV 2 and 
taking p = aip*ip; where a is a mass parameter, different from b, V = p/m and 
p is the momentum. 

Using the relations from the Hamilton- Jacobi theory 

± = e **<*)/* p = VS(x)=mV, (3) 
we can express the square of the velocity in terms of ip and ip* as follows, 

V = -i— Vln^, (4) 
so the energy-density becomes 



1 , ^ ,o ah 2 „ ^ ib ^ ib* 

oH^ 2 = ^^*Vln^-Vln^-, 5 

1 Sm z ib* ib 
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from which we immediately conclude that the corresponding nonlinear potential 
term associated with the kinematical pressure term is 

ah 2 - ib - ib* 

— ^Vln-f- • Vln-V. 6 

8m 2 tp* ip 

Hence a candidate for a NLSE is 

iB f " + u * ~ Wr*M + 5? ■ f '* T ) *■ (7) 

Here the Hamiltonian is Hermitian and a ^= b both are mass-energy parameters 
to be determined experimentally. As far as we know, this NLSE has not been 
derived so far. 

The new term can be written also in the form 

Vln^.Vln^ = -fvin^y. (8) 

For simplicity purposes, from now on we shall refer to these nonlinear po- 
tential kinematic pressure terms as those terms of the type 

U km ~ (Vln^) 2 . (9) 

The reason we choose to impose that notation will become clear in the next 
sections. 

Our goal now is to derive NLSE directly from the fractal space time dynamics 
of a particle undergoing a Brownian random motion. And such fractal space 
time interpretation does not require to add ad-hoc terms like 

1. David Bohm's quantum potential |^] into the Hamilton equation. See 
also §. 

2. To use the hydro-dynamical models discussed so far [|[ ||, [l(| nor to 
add the hydrostatic pressure and kinematic pressure terms to the Euler-Newton 
equations of motion as we have shown above, in an ad-hoc fashion. 

The new NLSE can be obtained from first principles if, an only is, we as- 
sume a fractal trajectory of a point particle associated with a Brownian random 
motion Q|. 

Before we begin, we deem it very important to add some more comments 
about the kinematic pressure terms 

V 2 ^-|Vln^| 2 , (10) 
2 2m m 

versus the hydrostatic pressure terms 

[— & -61n(V*^) (11) 
J P 

in the new NLSE. 
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The hydrostatic term —bln{^*rp)^ explicitly breaks homogeneity tj) — ► Xtp 
of the NLSE. Whereas the kinematic pressure term (U 2 ) / (2m){a/m)\V In^ip 
does preserve the homogeneity condition and the NLSE should scale with a A 
factor, fact that can be easily verified. 

The hydrostatic pressure term is not compatible with the motion kinematics 
of a particle executing a fractal Brownian motion. Only in the m — > oo limit 
(heavy particle) it is reasonable to speak of the static limit. See for a discus- 
sion of this limit. There are another deeper reasons to derive the NLSE from an 
underlying dynamics of a particle in a fractal space time, or more simply, from 
a fractal trajectory in a fixed space time background. We may, or may not, be 
switching on the quantum gravitational aspects of space time, so deeply linked 
to the non-linearity of QM. Perhaps the nonlinearity of QM is deeply inter- 
wined with the quantum gravitational aspects of a Cantorian-fractal-space-time 
§& 

Our goal is far less ambitious. Returning to our main points, the two NLSE 
have introduced two additional parameters of mass-dimension a, b (or one pa- 
rameter in the special case a = b). Such parameters need to be found exper- 
imentally. The advantage of the fractal formalism is that we will be able to 
relate the a, b parameters to the Planck constant ft itself, rather to have new 
parameters in physics unrelated to K. 



3 QM as mechanics in non different iable spaces 

We will be following very closely Nottale's derivation of the ordinary Scrodinger 
equation The readers familiar with this work may omit this section. Re- 
cently Nottale and Celerier Q following similar methods were able to derive the 
Dirac equation using bi-quaternions and after breaking the parity symmetry 
dx^ <-> —dx* 1 , see references for details. Also see the Ord's paper Q and the 
Adlers's book on quaternionic QM For simplicity the one-particle case is 
investigated, but the derivation can be extended to many-particle systems. 

In this approach particles do not follow smooth trajectories but fractal ones, 
that can be described by a continuous but non-differcntiable fractal function 
f(t). The time variable is divided into infinitesimal intervals dt which can be 
taken as a given scale of the resolution. If $(t,t',dt) is a smoothing function 
centered on t, for example a step function of width 2dt, a continuous and dif- 
ferentiable approximation to the true fractal r(t) can be constructed as follows, 

oo 

f(t,dt)= J $(t,t',dt)r(t')dt'. (12) 

— OG 

While r(t) = r(t, 0) is non-differentiable, any r(t, dt), called "fractal trajectory" , 
is differentiable for all dt ^ 0. 

Non-differentiability implies a lost of causality. For this reason the fractal 
trajectories are good candidates to describe the quantum behavior. Feynman in 
his path integral formulation of QM already found a interesting result involving 
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the time scale dt: When seen at a time scale dt, the quantum mechanical mean 
quadratic velocity of a particle is (v 2 ) oc 1/dt. This can be easily explained 
by an argument involving fractals. If the trajectory is a fractal curve of fractal 
dimension D, the space and time resolutions are related by dt oc dx D , so that 
(v 2 ) oc (dx/dt) 2 oc dt 2 ^ 1 ' ' D ^ 1 \ The comparison with Feynman's result leads to 
D = 2. 

Here we omit the details of the arguments leading to dt oc dx D , which the 
interested reader can find at jlj and references therein. They has as ingredients 
the scale dependence of any fractal curve and the renormalization group. 

A fractal function f(x,e) can have, besides the derivative df(x,e)/dx, the 
new derivative with respect to the scale, df (x, e)/de. It was found useful to use 
In e instead of e as the variable for resolution. Renormalization group arguments 
say that the following relation is valid fj| , 

^A=a(x)+bf(x,e), (13) 
a me 

this means that the variation of / under an infinitesimal scale transformation 
dine depends only on / itself. This differential equation can be integrated to 
give us 



f(x,e) = f (x) 



l + ((x) - 



(14) 



X~ b ((x) is an integration constant and fo(x) = ~a(x)/b. This result says that 
any fractal function can be approximated by the sum of two terms, one inde- 
pendent of the resolution and other resolution dependent. Due to the resolution 
dependence is associated to the fractal properties, and those are product of the 
non-differentiability, then is expected that Ci x ) is a fluctuating function with 
zero mean. 

Provided than a ^ and b < two cases can be considered: (i) e <C A, the 
scale dependent term is dominant and / is given by a scale-invariant fractal-like 
power law with fractal dimension D = b — 1, namely f(x,e) — fo(x)(\/e)~ b . 
(ii) if e ^> A then / becomes independent of the scale. A is the de Broglie wave 
length. 

A continuous but non-differentiable function f(t) at t has two possible values 
of the derivative at t, for this reason its approximation by a fractal function 
requires considering "left" and "right" derivatives. For the position vector the 
following two infinitesimal differences can be considered, 



fit + dt,dt) - r(t,dt) = b+ (r, t)dt + £+ (t , dt) — 



dt^ 



f(t, dt) - r(t - dt, dt) = 6_ (r, t)dt + £_ (t, dt) I — ] 



T oJ 



(15) 



where j3 — 1/D, and b + and 6_ are average forward and backward velocities Q. 
Adopting the non standard analysis formulation, dt is also the time scale. 
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The instantaneous velocities are easily obtained from equation (filf), 

v ± (f,t,dt) = b ± (f,t)+^ ± (t,dt)l-\ , (16) 

In the quantum case, D = 2, then [3 = 1/2, so that dt^~ x is a divergent quantity, 
from which is evident the non-differentiability. 

Then, following the definitions given by Nelson in his stochastic QM ap- 
proach (Lemos in Jl6| p. 615; see also jL7|, |ll|), Nottale define mean backward 
an forward derivatives as follows, 

^±M= hm (r(t + ^-m\ 
dt At^±o \ At I 

from which the forward and backward mean velocities are obtained, 



(18) 



s 


■I 


fd+ 


d\ 


.1 (d+ 


d_ 


dt 




{ dt 


"dt) 


l 2 \ dt 


dt 



For his deduction of Schrodinger equation from this fractal space-time clas- 
sical mechanics, Nottale starts by defining the complex-time derivative operator 

(19) 

which after some straightforward definitions and transformations takes the fol- 
lowing form, 

— = — + V -V -iDV 2 . (20) 
dt dt y ' 

D is a real- valued diffusion constant to be related to the Planck constant. Now 

we are changing the meaning of D, since no longer a symbol for the fractal 

dimension is needed, it will have the value 2. 

The D comes from considering that the scale dependent part of the velocity 

is a Gaussian stochastic variable with zero mean, (see de la Pefia at (l6) p. 428) 

(d&id&y) = ±2DS tJ dt. (21) 

In other words, the fractal part of the velocity £, proportional to the C, amount 
to a Wiener process when the fractal dimension is 2. 

Afterwards, Nottale defines a set of complex quantities which are generaliza- 
tion of well known classical quantities (Lagrange action, velocity, momentum, 
etc), in order to be coherent with the introduction of the complex-time deriva- 
tive operator. 

The complex time dependent wave function ip is expressed in terms of a 
Lagrange action S by tp — e lS /( 2mD ). S is a complex- valued action but D is 
real-valued. The velocity is related to the momentum, which can be expressed 
as the gradient of S, p = VS. Then the following known relation is found, 

V = —2iDVlnip. (22) 
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The Schrodinger equation is obtained from the Newton's equation (force = 
mass times acceleration) by using the expression of V in terms of the wave 
function t/j, 

S 8 
- W = m— V = -2imD — Vlntb. (23) 
dt dt Y v ' 



Replacing the complex-time derivation (|20|) in the Newton's equation gives 

us 

- W= -2im f-D-^Vln^J -2DVyD^-j . (24) 

Simple identities involving the V operator were used by Nottale. Integrating 
this equation with respect to the position variables finally yields 

° 2 VV + i Df -iL« = 0, (25) 

up to an arbitrary phase factor which may set to zero. Now replacing D by 
h/(2m), we get the Schrodinger equation, 

The Hamiltonian operator is Hermitian, this equation is linear and clearly is 
homogeneous of degree one under the substitution ip — > Xip. 



4 Nonlinear QM as a fractal Brownian motion 
with a complex diffusion constant 

Having reviewed Nottale's work |Q we can generalize it by relaxing the assump- 
tion that the diffusion constant is real; we will be working with a complex-valued 
diffusion constant; i.e. with a complex-valued h. This is our new contribu- 
tion. The reader may be immediately biased against such approach because the 
Hamiltonian ceases to be Hermitian and the energy becomes complex- valued. 
However this is not always the case. We will explicitly find plane wave solutions 
and soliton solutions to the nonlinear and non-Hcrmitian wave equations with 
real energies and momenta. 

For a detailed discussion on complex-valued spectral representations in the 
formulation of quantum chaos and time-symmetry breaking see Also a 

complex-valued time and two-times (see |22j] and references therein) complex- 
valued dimensions have been discussed in ]20| , [2l| | . 

Nottale's derivation of the Schrodinger equation in the previous section re- 
quired a complex-valued action S stemming from the complex- valued velocities 
due to the breakdown of symmetry between the forwards and backwards ve- 
locities in the fractal zigzagging. If the action S was complex then it is not 
farfetched to have a complex diffusion constant and consequently a complex- 
valued h (with same units as the complex- valued action). 
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Our derivation follows closely that of Nottale, sketched into the previous 
section, but with some crucial differences in the evaluation of the correlation 
functions and the definition of the complex-time derivative operator, respec- 
tively. 

Before the derivation further comments on complex-energies are in order. 
The energy functional Eqm contains imaginary components. Since meaningful 
physical solutions demand real-valued energies this imposes constraints on the 
physically acceptable states in these non-linear QM equations, see Puszkarz |pd| . 

Complex energy is not alien in ordinary linear QM. They appear in optical 
potentials (complex) usually invoked to model the absorption in scattering pro- 
cesses Q and decay of unstable particles. Complex potentials have also been 
used to describe decoherence Q . The accepted way to describe resonant states 
in atomic and molecular physics is based on the complex scaling approach, which 
in a natural way deals with complex energies |^4| . We will show that real- valued 
energy solutions exist to the NLSE based on a fractal Brownian motion. 

The imaginary part of the linear Schrodinger equation yields the continuity 
equation dp/dt + V • (pV) = 0. Because, as we shall see, our potential is 
complex, the imaginary part of such potential acts as a source term in the 
continuity equation. 



Before, Nottale wrote, 



(d(±d( ± ) = ±2Ddt, (27) 



with D and 2mD — h real. 

Now we set 



(d(±d(±) = ±{D + D*)dt, (28) 

with D and 2mD — Ti — a + i(3 complex. 

The complex-time derivative operator becomes now 

* = ! + ?. V-i( D + D.)V*. (29) 

In the real case D = D*. It reduces to the complex-time-derivative operator 
described previously by Nottale. 

Writing again the ip in terms of the complex action S, 

i> = e *S/(2mD) = jS/h^ (3Q) 

where S, D and h are complex-valued, the complex velocity is obtained from 
the complex momentum p — VS as 

V = -2iDVln^>. (31) 

The NLSE is obtained after we use the generalized Newton's equation (force 
= mass times acceleration) in terms of the ip variable, 

S 8 
— V£7 = m — V — ~2imD — Vlnib. (32) 
dt dt J 
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Replacing the complex-time derivation ( |2^ ) in the generalized Newton's 
equation gives us 



VU = 2im 



D^Vlnip - 2iD 2 (yiniP ■ V)(Vln^) - ^(D -!- LV )DY 1 [\\u r) 



(33) 

Now, using the three identities (1) VV 2 = V 2 V, (2) 2(Vln^ • V)(Vln^) = 
V(Vln-0) 2 and (3) V 2 ln-0 = V 2 -0/V> - (Vln^) 2 allows us to integrate such 
equation above yielding, after some straightforward algebra, the new NLSE 
that has the nonlinear (kinematic pressure) potential found before (pi), 



Note the crucial minus sign in front of the kinematic pressure term and that 
h = a + id = 2mD is complex. When 8 = we recover the linear Schrodinger 
equation. 

The nonlinear potential is now complex-valued in general. Defining 

^="^M . ( 35 ) 



2m fi 

and U the ordinary potential, then the NLSE can be rewritten as 



^dtp ( U 2 a 2 



in it={-^ +u+m r- (36) 

This is the fundamental nonlinear wave equation of this work. It has the form 
of the ordinary Schrodinger equation with the complex potential U + iW and 
the complex fi. The Hamiltonian is no longer Hermitian and the potential it- 
self depends on %p. Nevertheless one could have meaningful physical solutions 
with real valued energies and momenta. Like the plane-wave and soliton solu- 
tions. Notice that the new NLSE obeys the homogeneity condition -0 — > At[> 
for any constant A. All the terms in the NLSE are scaled respectively by a 
factor A. We did not obtain the hydrostatic pressure term — 6(ln'0*V')' ! /' which 
breaks the homogeneity condition for a simple reason: We are studying the true 
kinematics and dynamics of a particle of mass m undergoing a fractal Brow- 
nian motion. It would be meaningless to have a hydrostatic pressure term in 
such a model. Moreover, our two parameters a, 8 are intrinsically connected 
to a complex Planck constant Ti = a + id rather that being ah-hoc constants 
to be determined experimentally. Thus, the nonlinear QM equation derived 
from the fractal Brownian motion with complex-valued diffusion coefficient is 
intrinsically tied up with a non-Hermitian Hamiltonian and with complex- valued 
energy spectra [|l3|. To be more precise, the nonlinear 8 term in ([36]) is really 
the nonlinear partner of the kinetic energy term. 

We will show that despite having a non-Hermitian Hamiltonian wc still could 
have eigenfunctions with real valued energies and momenta. When h is real 
(8 = 0) and the NLSE is linearized back to the ordinary one. 
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The reader may ask why not simply propose as a valid NLSE the following, 



ih l£ = + W+ ^--|Vln^|V. (37) 

ot 2m 2m m 

Such equation (with the ordinary real Planck constant) is based on a real Hamil- 
tonian that satisfies the homogeneity condition. It also admits soliton solutions 
of the type 

^ = CA[x-Vt)e i{ - kx - ut \ (38) 

A(x — Vt) is a function to be determined by solving the NLSE. Galilean invari- 
ance imposes that the soliton is a traveling wave, a function of x — Vt. We will 
present explicit expressions for the function A(x — Vt) afterwards. Therefore, 
in principle, this NLSE is a more suitable candidate that the Bialynicki-Birula 
and Mycielski NLSE with a nonlinear potential term —b\n(ip*ip) that breaks 
homogeneity and introduces also a tj)* dependence into the wave equation. 

The only problem (perhaps there are others) with the NLSE above is that it 
suffers also from an extraneous dispersion relation. Plugging-in the plane-wave 
solution if) ~ e - % i Et ~P x )/ h one g e t s an extraneous energy-momentum relation, 
after setting U = 0, 

E = f(l + ^), (39) 
2m V to/ 

not the usual E = p e /(2m). So in this case we have that Eqm ^ Eft (FT 
means field theory). 

It has been known for some time, see Puskarz |n]], that the expression for 
the energy functional in nonlinear QM does not coincide with the QM energy 
functional, nor it is unique. The simplest way to see this is, for example, writing 
down the Birula and Mycielski NLSE (^) in the Weinberg form Q 

at ~ or ' 1 j 

where ip an d "0* are a pair of canonically-conjugate variables. The real-valued 
Hamiltonian density is given by 

h 2 

HN), ib*) = -—i)*V 2 ib + Ui>*i) - 6V* H^ibU + bib*ip, (41) 
2m 

using Eft = J d 3 rH, so we can see it is different from (H)qm- Notice the last 
term in @. Hence, one can immediately sec that the H(ip, ib*) stemming from 
the field theory approach does not coincide with the Birula-Mycielski Hamil- 
tonian. They differ by a constant, Eft — Eqm = J d 3 rbib*ip = b. Exactly 
like it occurs when we plug in the plane wave solution into the NLSE with the 
nonlinear potential, real valued, 
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Notice that this problem does not occur in the fractal-based NLSE, because such 
NLSE is written entirely in terms of the i/j variable and does not contain the t/j* 
variable explicit or implicitly, like it occurs in the Birula-Mycielski NLSE. 

The classic Gross-Pitaveskii NLSE (of the I960'), based on a quartic in- 
teraction potential energy, relevant to Bose-Einstein condensation, contains the 
nonlinear cubic terms in the Schrodinger equation, after differentiation, (ip*ip)ip. 
This equation does not satisfy the Weinberg homogeneity condition and also the 
Eft differs from the Eqm by factors of two. 

In the fractal-based NLSE there is no discrepancy between the quantum- 
mechanical energy functional and the field theory energy functional. Both are 
given by 

*#£55 = "I^Vty + W " ^(VlnVOV (43) 

The NLSE is then unambiguously given by equation (|4(]), H(ip,ip*) is homoge- 
neous of degree 1 in A respect to ip. This is why we push forward the NLSE 
derived from the fractal Brownian motion with a complex-valued diffusion co- 
efficient. Such equation does admit plane-wave solutions with the dispersion 
relation E = p 1 /(2m). It is not hard to see that after inserting the plane wave 
solution into the fractal-based NLSE we get (after setting U = 0), 

2m h h 2 1 ft 2m 2m fr 2m ' ^ ' 

since h = a + i(3. So the plane-wave is a solution to the fractal-based NLSE 
(when U — 0) with a real-valued energy and which has the correct energy- 
momentum dispersion relation. 



5 Soliton solutions to the fractal based NLSE. 
One dimensional case 

Let us find soliton solutions to the fractal-based NLSE given by (HJ) , in the free 
particle case. We set the ansatz (one-dimensional for simplicity) 

i/f = CA(x - Vt)e~ {Et ~ px)/h . (45) 

The function A must be complex-valued, otherwise no real-valued energy solu- 
tions exist. Then we set 

A{x-Vt)=F(x-Vt)+iG{x-Vt), (46) 

and plugging-in if) with this A into the fractal-based NLSE (J|J) yields 2 coupled 
differential equations, after separating the real and imaginary parts, respectively, 
which yield, in principle, the functions F and G. 

For example, the soliton solution to the NLSE with the —bln(ip*tp) is of the 
form |, 

$(x, t) = Ce a / B e-( B /W x - Vt+ V 2 e ikx ~ iut , (47) 
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where c, a, B, d are numerical constants which can depend on h, m and b. 

As mentioned before, plane wave solutions to the NLSE based on the —b \n.{ij)*tp) 
potential exist but they have extraneous dispersion relations. For example, the 
energy-momentum relation turns out to be 0j E = Twj — p 2 /(2m) + 61n(27r). 
Thus, plane-wave solutions do not seem to have physically meaningful interpre- 
tation. This was another reason why we believe that this NLSE has problems. 
Besides, we remarked already that this NLSE breaks the homogeneity condition 
as well. 

To finalize we will find the soliton solutions to the NLSE based on the kine- 
matic pressure potential | V In | 2 terms, given by equation (j37|), in the free 
particle case U = 0. 

Earlier on we have shown that it admits plane wave solutions with the ex- 
traneous dispersion relation E = p 2 /(2m)(l + a/m). It obeys the homogeneity 
condition: Under scaling of i/j by A the NLSE scales with an overall factor of A 
as expected. 

Notice that if we wish to have a Hermitian Hamiltonian we must take the 
absolute value | V In ^ | 2 instead of (Vln^) 2 for our potential. Notice this impor- 
tant difference between these non-linear potentials in the fractal-based NLSE 
versus the kinematic pressure based one. 

Pluggin-in the ansatz ip = CF(x — Vt)e~^ Et ~ px ^ n into the kinematic pres- 
sure NLSE for the free particle case U — yield for the imaginary parts 
—ihVF' = —i(Ti/m)F'p. Then, for any F we have V — p/m. Therefore the 
ansatz is consistent with the de Broglie relations p = hk and E — hto, as ex- 
pected from this NLSE soliton solution. 

The real parts give the differential equation 



a ,„, s0 1 



a 



F"F (F') 2 + - s 2mE -p 2 {l + — ) F = 0. (48) 



in 



The solutions to this nonlinear differential equation yield F(x — Vt). This 
differential equation involves the derivatives F' , F" and is much harder to solve 
than the differential equation given in Q that involves F" but not F' . 

For example a nonlinear differential equation which involves F" but not F' 
is F" — F 3 — 0. Such equation, after multiplying both sides by F', can then be 
integrated by quadratures, J dF/F 2 = j dy/2 1 / 2 . 



6 Concluding remarks 

Based on Nottale and Ord's formulation of QM from first principles; i.e. from 
the fractal Brownian motion of a massive particle we have derived explicitly a 
nonlinear Schrodinger equation. Despite the fact that the Hamiltonian is not 
Hermitian real-valued energy solution exist like the plane wave and soliton so- 
lutions in the free particle case. The hydro-dynamical model analog of this 
fractal-based NLSE yields another new NLSE with Hermitian (real) Hamilto- 
nian. The remarkable feature of the fractal approach versus all the nonlinear 
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QM equation considered so far is that the quantum mechanical energy functional 
coincides precisely with the field theory one. 

The hydro-dynamical-based NLSE has a nonlinear (real) potential term 

ah 2 -> , tb -± , ib* , , . 

8^ V1 V' Vln ^ ™ 

with a the mass-energy parameter, bears a very rough similarity to the Starusz- 
kiewicz imaginary potential term in three dimensions 

- l-^r- < 50 > 

see pH , this potential is imaginary, 7 is a constant. 

The fractal model based NLSE admits plane wave (soliton solutions also) 
with the correct dispersion relation E — p z /(2m), real. Soliton solutions, with 
real-valued energy (momentum) are of the form 

ib ~ [F(x - Vt) + iG(x - V t)]e ipx/h - lEt/h , (51) 

with F, G two functions of the argument x — Vt obeying a coupled set of two 
nonlinear differential equations. 

It would be interesting to study solutions when one turns-on an external 
potential U 7^ 0. 

The reader may ask why concentrate on a complex diffusion constant to 
generate a nonlinear Schrodinger equation with a non-Hermitian Hamiltonian 
when one could have written from the start a NLSE with a Hermitian Hamilto- 
nian that obeys the Weinberg homogeneity conditions and also with the correct 
energy dispersion relations. 

Starting from the fundamental equations: ib = e lS ^ So = e lS ^ h ° and the 
generalized Newtonian law, written in terms of the Nottale complex derivative 
operator and the ip: 

W = iS [ dV !, n ^ ~ — (Vln^ ■ V)(VlnV) + £>V 2 (V In </;)}], (52) 
at m 

after adding and subtracting the quantity Z?oV 2 (Vln-0), and using the 3 vector- 
calculus identities used in (|36|), we get a nonlinear correction to the Schrodinger 
equation 

^(h-h a )(V 2 \n^, (53) 
Zm 

where So = h = 2mD and h = 2mD 7^ 2rnZ? - 

As desirable as this NLSE may look a close inspection reveals that the non- 
linearity is just an artifact of the definition of ib. It looks nonlinear from the ib 
perspective. It is not difficult to see that under a re-definition of the wavefunc- 
tion 

if,' = e zS/h = e * s /( 2mD ) (54) 
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the generalized Newtonian law of motion of a particle undergoing a fractal 
Brownian motion given by ( [52] ) yields now the standard linear Schrodinger 
equation for the ip' wavefunction. Where we have chosen for a new value of 
S' = h = 2mD. 

It is important to emphasize that the diffusion constant is always chosen to 
be related to Planck constant as follows: 2mD = Ti which is just the transition 
length from a fractal to a scale-independence non-fractal regime discussed by 
Nottalc in numerous occasions. In the relativistic scale it is the Compton wave- 
length of the particle (say an electron): A c — h/ (mc). In the nonrelativistic case 
it is the de Broglie wavelength of the electron. 

Therefore, the NLSE based on a fractal Brownian motion with a complex 
valued diffusion constant 2mD = h = a + i/3 represents truly a new physical 
phenomenon in so far as the small imaginary correction to the Planck constant 
(unobserved in present day experiments) is the hallmark of nonlinearity in QM. 
For other generalizations of QM see experimental tests of quaternionic QM 
(in the book by Adler pgj] ). Equation ( j36| ) is the fundamental NLSE of this 
work, where the /3 term is essentially the nonlinear partner of the linear kinetic 
energy term in comparison to all other approaches which focused on nonlinear 
modifications of the potential. 
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